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Abstract 

A search for first generation scalar and vector leptoquarks (LQ) as well as for 
squarks (q) in R-parity violating SUSY models has been performed using e + e~ 
collisions collected with the OPAL detector at LEP at an e + e~ centre-of-mass en- 
ergy t/s^T of 189 GeV. The data correspond to an integrated luminosity of about 
160 pb _1 . The dominant process for this search is eq — ► LQ/q — ► eq, uq, where a 
photon, which has been radiated by one of the beam electrons, serves as a source 
of quarks. The numbers of selected events found in the two decay channels are 
in agreement with the expectations from Standard Model processes. This result 
allows to set lower limits at the 95 % confidence level on the mass of first gener- 
ation scalar and vector leptoquarks, and of squarks in R-parity violating SUSY 
models. For Yukawa couplings A to fermions larger than y/4iia em , the mass limits 
range from 121 GeV/c 2 to 175 GeV/c 2 (149 GeV/c 2 to 188 GeV/c 2 ) depending 
on the branching ratio (3 of the scalar (vector) leptoquark state. Furthermore, 
limits are set on the Yukawa couplings A for leptoquarks and A' ljfc for squarks, 
and on (3 as a function of the scalar leptoquark/squark mass. 
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1 Introduction 



Leptoquarks (LQ) are coloured spin or spin 1 particles carrying both baryon 
(B) and lepton (L) quantum numbers. They appear in many extensions of the 
Standard Model as a consequence of the symmetry between the lepton and quark 
sectors. The Buchmiiller-Ruckl-Wyler (BRW) model Q used in this paper as- 
sumes lepton and baryon number conservation. Moreover the simplifying as- 
sumption is made that a given leptoquark couples to just one family of fermions 
which means that only first generation leptoquarks can be produced in electron- 
photon scattering. The first generation leptoquarks may decay into either an 
electron^ and a quark or into a neutrino and a quark. The branching ratio of the 
decay into an electron and a quark is commonly denoted by (3. Table |l] shows all 
the possible leptoquark states considered in this paper (scalar and vector) using 
the most common nomenclature || along with their electric charge and fermion 
number F = L + 3B. The branching ratios (3 given in this table assume that 
the leptoquarks couple only to Standard Model particles. To respect the existing 
limits on the product Al • \r we assume that for the scalar states S and Si/ 2 and 
the vector states V and Vi/ 2 either the left-handed or the right-handed coupling 
must vanish, i.e. Al • Ar = 0. 

At an e + e~ collider different diagrams are expected to contribute to single 
leptoquark production in electron-photon collisions j|, || || [7|, ||. The dominant 
diagram is eq — > LQ (Fig. [I]), where a photon, which has been radiated by one 
of the beam electrons, serves as a source of quarks through its fluctuations into 
hadronic states. The electron-quark interaction produces a leptoquark which is 
assumed to decay subsequently into an electron or a neutrino, and a quark. The 
photon remnant may disappear down the beam-pipe or add some activity in the 
forward region of the detector. The diagrams shown in Fig. |^ are also relevant, 
whereas e + e~ annihilation diagrams with a single leptoquark radiated in the final 
state and diagrams with quark and/or leptoquark exchange in the t-channel but 
without photon exchange are suppressed. The signature of single leptoquark 
events is one hadronic jet with high transverse momentum, balanced either by 
one isolated electron or by missing transverse energy due to the neutrino. Both 
topologies are studied in this paper. 

Squarks (q) in supersymmetric models with R-parity violation have the same 
production mechanism as some leptoquarks. R-parity is a quantum number which 
equals +1 for particles and —1 for their superpartners. Table |I] shows the cor- 
respondence between the squark and the leptoquark states. R-parity conserving 
decays are possible for squarks, in addition to the R-parity violating leptoquark 

1 Charge conjugation is implied throughout this paper for all particles, e.g. positrons are also 
referred to as electrons. 
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decay modes. The ratio between the R-parity conserving and violating modes 
depends on the parameters of supersymmetry and on the size of the coupling. For 
this analysis the branching ratio for R-parity conserving decays has been set to 
zero and consequently (3 has the same value for squarks and the corresponding lep- 
toquarks. Supersymmetry allows only left-handed couplings to leptons for these 
states. In the most general superpotential of the Minimal Supersymmetric Model 
(MSSM), the renormalizable gauge invariant operator which describes the cou- 
pling of squarks to quarks and leptons and violates R-parity is X'^L^Q-^D^ ||, 
where i, j and k are generation indices of the left-handed doublet superfields 
of leptons (Ll) and quarks (Qi), and right-handed singlets of down-type quarks 
(_D R ). Only couplings A' ljfc to first generation leptons are considered in this paper. 

Several experiments have searched for leptoquarks. DELPHI analysed the 
single scalar and vector leptoquark production setting limits on the mass ranging 
from 134 GeV/c 2 to 171 GeV/c 2 at an e + e~ centre-of-mass energy \/s ee of 183 
GeV [[| and for a Yukawa coupling A > ^Ana cm . The HI collaboration at 
HERA has searched for leptoquarks in deep-inelastic neutral current (NC) and 
charged current (CC) electron-proton scattering at high Q 2 ||10|| . First generation 



scalar (vector) leptoquarks have been excluded for masses up to 275 GeV/c 2 
(284 GeV/c 2 ) for A > ^4na em and fermion number F — 0, while masses up to 
about 200 GeV/c 2 have been excluded for \F\ = 2 leptoquarks. Leptoquark limits 
as a function of the assumed couplings have also been obtained from fermion pair 
production at LEP2 at ^ ee = 130 - 183 GeV 



Leptoquark pair production limits have been obtained by the LEP experi- 
ments at v^ ee = M z and at ^/s ee = 183 GeV @, and by the DO and 



CDF experiments JL5[] for leptoquarks of the three generations. The leptoquark 
pair production limits are independent of the Yukawa coupling A in pp scat- 
tering. In e + e~ scattering the leptoquark pair production cross-section can be 
considered independent of A only in the region of small A where t-channel quark 
exchange can be neglected compared to the s-channel diagram. At LEP2 energies, 
\/s ee = 183 GeV, the mass limits for first generation leptoquarks vary between 
80 GeV/c 2 and 90 GeV/c 2 , depending on the leptoquark state fl3fl . Scalar lepto- 
quarks of charge —1/3 are only excluded for Mlq < Mz/2 |T3]. A combination 
of the CDF and DO search results for a first generation scalar leptoquark yields 
a lower mass limit of 242 GeV/c 2 for = 1 [Pi The CDF and DO collabora- 
tions [I7J as well as the four LEP experiments [18| have also searched for pair 



production of R-parity violating squarks. 

Leptoquark pair production limits obtained at LEP are sensitive to the mass 
region Mlq < A /s^/2, whereas single leptoquarks can be produced almost up 
to the kinematic limit, Mlq = y/s^. Even though the leptoquark mass range 
covered by CDF and DO in pp scattering is higher for most leptoquark states, the 
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Table 1: All possible scalar (S) leptoquarks/squarks and vector (V) leptoquarks 
in the BRW model with their electric charge in units of e, their fermion number 
F, their decay modes and the corresponding branching ratio f3 for the decay into 
an electron and a quark. 



6 



analysis presented here is more sensitive in the low (3 region {(3 — > 0). For (3 = 
no production in eq collisions is possible. In addition, this search is also sensitive 
to the production of ec (and es) states, but only flavour-diagonal couplings A are 
considered in this paper. 

We present a search for leptoquarks with M L q > 80 GeV in electron-photon 
scattering using data corresponding to an integrated luminosity of 164.7 pb _1 
(eq channel) and 158.4 pb _1 (z/q channel) at e + e~ centre-of-mass energies of 189 
GeV. 



2 The OPAL detector 



The OPAL detector is described in detail in |19| . It is a multipurpose appara- 
tus having nearly complete solid angle coverage with excellent hermeticity. The 
central detector consists of two layers of silicon micro-strip detectors f2(| sur- 



rounding the beam-pipe and a system of gaseous tracking chambers inside a 
0.435 T solenoidal magnetic field. 

The lead-glass electromagnetic calorimeter (ECAL) with a presampler is lo- 
cated outside the magnet coil. It provides, in combination with the forward 
calorimeters (FD), the forward scintillating tile counter (the "MIP plug") [ pl|j , 
and the silicon-tungsten luminometer (SW) ||22|| , a geometrical acceptance down 
to 25 mrad from the beam direction. The SW luminometer measures the inte- 
grated luminosity using small-angle Bhabha scattering events |23fl. The magnet 
return yoke is instrumented for hadron calorimetry (HCAL). It is surrounded by 
several layers of muon chambers. 



3 Kinematics and Monte Carlo simulations 



The Monte Carlo simulation of the process e + e — > LQ + X is done with the 
program ERATO-LQ which can generate all the states listed in Table [I] and 



calculates the cross-sections for the scalar and vector states^. 

The total cross-section for the production of leptoquarks of mass Mlq can be 
written as a convolution of the probability to find a photon with the momentum 
fraction z in the electron, approximated here by the Weizsacker- Williams effec- 
tive photon distribution f^/^z) []2"4| , and the probability to find a quark in the 



2 All total cross-sections in this paper are defined as a sum of the particle and the anti-particle 
state. 
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photon. This probability can be parametrised by parton distribution functions 
(pdf) fq/ y (x, /i 2 ) of the photon, which are evaluated at the scale /i 2 = M 2 ||. 
The Bjorken scaling variable x is given by x = M 2 /zs. With these assumptions 
the total cross-section for scalar leptoquark production is: 



a(e + e - _> L Q + X) = ^ f -f, /e (z)f qh (M 2 /(zs), M 2 ) (1) 

where f q /.y(M 2 /(zs), M 2 ) is obtained by convoluting the parton level cross sec- 
tion with the quark distribution / q / 7 (z, M 2 ) in the photon. In case of unpolarised 
electron beams, the total cross-section for the production of vector leptoquarks is 
twice as large as the cross-section for scalar leptoquarks |JJ . This approach based 
on the pdf is used in the calculations of the diagram shown in Fig. [I] by Donch- 
esky PJ and it is also implemented in the Monte Carlo generator PYTHIA 



However, PYTHIA can only generate scalar leptoquarks with charge —1/3. It 
has therefore only been used to check the cross-sections calculated by ERATO- 
LQ p5 |. This Monte Carlo generator uses a perturbative calculation of the dia- 



grams in Fig. m and 0. It is expected to give the correct angular distributions 
both for the scalar and vector leptoquarks. 

Cross-sections calculated with ERATO-LQ and PYTHIA using A = v / 4vra cm 
and charge —1/3 (e.g for So) are shown in Table Q Using a different parametri- 
sation of the pdf, GRV [fZ7[] instead of SaS-lD p8[, has almost no effect on the 
resulting cross-section. In ERATO-LQ the exact total cross-section given in 0, 
taking into account all possible diagrams, can also be calculated. In comparison, 
the total cross-section calculated perturbatively by the MC generator ERATO- 
LQ is 10 — 20% smaller than the exact total cross-section used to calculate the 
limits. 

ERATO-LQ generates the four-vectors of the direct decay products of the 
leptoquark. A scalar leptoquark decays isotropically in its rest frame leading to 
a flat distribution in the variable y — (1 + cos#*), where 6* is the decay angle of 
the lepton relative to the incident quark in the leptoquark centre-of-mass frame. 
The decay angles in the decays of vector leptoquarks are distributed according to 
da/dy oc (1 — y) 2 . The simulated photon is always real, i.e., the negative squared 
four- momentum of the photon, Q 2 , is identical to zero. 

The partial decay widths of a scalar (S) and a vector (V) leptoquark are 

3 A 2 

r s = -r v = — m lo . (2) 

b 2 v 16vr g K J 

Since the leptoquark carries colour, it could hadronise before its decay into 
fermions. This effect is taken into account in the systematic uncertainties but not 
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Table 2: The total cross-section for the single production of the state So 
with charge —1/3 as a function of the leptoquark mass Mlq calculated by the 
Monte Carlo generators ERATO-LQ and PYTHIA at y/s ee = 189 GeV using 

A = y/A7Ta em . 



in the standard Monte Carlo simulation, since it should only be important for 
r s ,v < Aqcd- The decay width T s is 16 MeV for A = 0.1 and M LQ = 80 GeV/c 2 . 



JETSET |26| is used to perform the hadronisation of the leptoquark decay 
products. It has been checked with ERATO-LQ that the event properties of the 
different leptoquark states are very similar which allows to simplify the generation 
considerably: For each of the seven masses listed in Table |^, samples of 3000 
events for scalar and vector leptoquarks were generated for the two leptoquark 
decay modes separately and only for one state. Also, no extra squark events 
needed to be generated. 

All relevant Standard Model background processes have been studied using 
Monte Carlo generators. Multi-hadronic events (e + e~ — > qq(7)) have been sim- 
ulated with PYTHIA 5.722 f2|]. KORALZ 4.02 J29] has been used to generate 
-> t + t~ and BHWIDE [RUJ to generate the Bhabha process 



the process e + e 
e + e~ — > e + e~. 



Deep inelastic e'j events in the range Q 2 > 4.5 GeV 2 including charged current 
deep inelastic scattering (CC DIS) events have been simulated with HERWIG 
5.8 f3lfl . PHOJET 1.10 |32] has been used to generate hadronic two-photon 
events (i.e. e + e~ — > e + e~ hadrons) in the range Q 2 < 4.5 GeV 2 . Leptonic two- 
photon events have been generated with Vermaseren [Q. The events generated 
with HERWIG, PHOJET and Vermaseren are called two-photon events. An 
alternative hadronic two-photon sample for systematic studies has been generated 
using F2GEN @ for Q 2 > 4.5 GeV 2 and PYTHIA for Q 2 < 4.5 GeV 2 . Other 
processes with four fermions in the final states, including W pair production, 
have been simulated with grc4f |J5 1 and an alternative sample has been generated 
with KORALW [S6|. The generated signal events and all the background Monte 
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Carlo events have been passed through a full detector simulation and the same 
reconstruction algorithms as the real data. 



4 Event analysis 



We search for events with one hadronic jet and either an electron or missing 
energy balancing the transverse momentum of this jet. The analysis uses tracks 
measured in the central tracking devices, clusters measured in the ECAL, the 
HCAL, the FD and the SW. In addition to quality requirements which ensure 
that the tracks have their origin close to the e + e~ interaction point, tracks must 
have more than 20 hits in the central jet chamber and more than half the number 
of hits expected for the given track. The transverse momentum of the track with 
respect to the beam direction must be greater than 120 MeV. Tracks with a 
momentum error larger than the momentum itself are rejected if they have fewer 
than 80 hits. Calorimeter clusters have to pass energy threshold cuts to suppress 
noise. To avoid double counting of particle momenta, a matching algorithm 
between tracks and clusters is applied p7 |. Clusters are rejected if the energy of 



the cluster is less than expected from the momentum of an associated track. If 
the cluster energy exceeds the expected energy by more than what is expected 
from the resolution, the expected energy is subtracted from the cluster energy. In 
this case the track momentum and the reduced energy of the cluster are counted 
separately. 

The tracks and remaining clusters are used as input to the jet finding al- 
gorithm and to determine the missing transverse energy of the event. Jets are 
reconstructed using a cone jet finding algorithm with a cone size R = 1 and a 
minimum transverse jet energy Et of 15 GeV |37]]. The cone size R is defined 



as R = y (Ar/) 2 + (A</>) 2 , with rj — — lntan(#/2) being the pseudorapidity, the 
azimuthal angle and 9 the polar angle in the laboratory frame in radians. At] 
and A(j) are the differences between the cone axis and the particle direction. No 
cut on the pseudorapidity of the jet is used at this stage. 



4.1 The electron plus hadronic jet channel 

In this search the selection cuts were optimised for a leptoquark decaying into a 
single quark and an isolated electron. The electron is identified by requiring a 
minimum of 20 hits used in the measurement of the specific energy loss, dE/dx, 
and a dE/dx probability for the electron hypothesis of more than 1%. Further- 
more, the ratio of the total energy of the electron measured in the ECAL to 
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the momentum of the track associated to this ECAL cluster must lie between 
0.7 and 2. The identified electron with the largest momentum was assumed to 
be the electron from the leptoquark decay. Since the electron is included in the 
jet search, it is usually reconstructed as a jet. Candidate leptoquark events are 
selected based on the following cuts, which are identical for scalar and vector 
lepto quarks: 

• The event must contain more than four tracks (n c h > 4). 

• Exactly two jets must have been reconstructed (nj et = 2). One of the jets 
must contain the highest energy electron. 

• The highest energy electron must have an energy E e greater than 2 GeV. 
The electron energy E e is the energy of the calorimeter clusters matched 
to the electron track. This cut is effective against all sources of Standard 
Model background, especially two-photon events. 

• The jet not containing the electron must consist of more than six particles 
(nqj > 6), where the number of particles is defined as the sum of the 
number of tracks and calorimeter clusters after matching. This cut reduces 
the number of e + e~ — > t + t~ and the number of Bhabha events and also 
some of the remaining leptonic two-photon events. 

• The total energy -Ehcal measured in the hadronic calorimeter has to be 
greater than 1 GeV. This cut is effective against Bhabha events. 

After this preselection, 5739 data events remain. The selection efficiencies are 
given in Table |^. They are significantly smaller for the vector leptoquark states 
than for the scalar leptoquark states due to the angular distribution of the de- 
cay electrons which is peaked at cos 8 = ±1 for vector leptoquark states. The 
following set of cuts is applied to the remaining events: 

EQ1) To ensure that most of the measured energy comes from the two jets, the 
energy of the hadronic jet and the electron energy E e must add up to 
more than 80% of the visible energy E vis , i.e. ((E^ + E c )/E vis > 0.8). This 
cut is efficient against multihadronic and four-fermion events (Fig. |3]a). 

EQ2) The ratio of the missing transverse energy ft? and the invariant mass My } 
of the two jets has to be fi^/M-^ < 0.15. The invariant mass My of the two 
jets is calculated from the four-vectors of the two reconstructed jets. This 
cut reduces the four-fermion background by a factor two (Fig. |]b). 
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EQ3) Since the electron is expected to be isolated, the difference between the 
electron energy E e and the energy E e j of the jet which contains the electron 
has to be smaller than 2 GeV (|E e — E e j| < 2 GeV). This eliminates most 
of the remaining e + e~ — > qq(7) events (Fig. |3]c). 

EQ4) The electron must lie in the angular region defined by | cos# e | < 0.8 (Fig. |3]d). 
This cut rejects mainly two-photon events with a scattered electron within 
the detector acceptance. 

In Table |3] the number of data events and the expected number of Standard Model 
background events taken from the Monte Carlo are shown after the preselection 
and after each subsequent cut. The number of Monte Carlo events has been 
normalised to the data luminosity. The selection efficiencies for three different 
scalar and vector leptoquark masses are also given. 



cuts 


Pre- 
selection 


F, ■ 

(EQl) 


(EQ2) 


E c — E ej 
(EQ3) 


cos 8 e 
(EQ4) 


M LQ [GeV/c 2 ] 


state 


80 


scalar 


49.6% 


44.7% 


39.4% 


37.9% 


28.4% 


120 


61.4% 


58.7% 


54.1% 


51.6% 


45.1% 


160 


70.7% 


70.5% 


65.4% 


60.2% 


55.0% 


80 


vector 


37.0% 


31.2% 


27.8% 


26.7% 


12.3% 


120 


50.4% 


47.5% 


43.9% 


41.6% 


29.4% 


160 


61.5% 


61.3% 


57.4% 


54.0% 


44.2% 


e + e~ — > qq 


4834 


200.1 


160.1 


2.7 


2.2 ±0.2 


e + e~ — ► t + t~ 


13.1 


5.2 


2.7 


0.7 


0.5 ±0.1 


e + e~ — > e + e~ 


0.5 


0.5 


0.5 


0.5 


0.2 ±0.1 


e + e~ — > 4 fermions 


817.2 


72.8 


33.0 


18.2 


12.8 ±0.5 


two-photon 


62.2 


30.0 


26.2 


25.5 


6.3 ±1.0 


total BG 


5727 


308.5 


222.4 


47.5 


21.9 ± 1.1 


data 


5739 


270 


194 


36 


21 



Table 3: Selection efficiencies for three different leptoquark masses for scalar and 
vector leptoquarks. The remaining number of data events and the expected num- 
ber of background (BG) events are also listed after each cut of the electron plus 
hadronic jet selection. The Monte Carlo background is normalised to the data 
luminosity of 164.7 pb _1 at 189 GeV. The errors on the Monte Carlo background 
are statistical. 

Figs. |^a-d show the distributions of some of the cut variables for data, Stan- 
dard Model background and the leptoquark state Vq with a mass of 120 GeV. 
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The data distributions are in general well described by the Monte Carlo simu- 
lation. After all cuts the Standard Model background is expected to be mainly 
due to four-fermion and two-photon interactions. 

The | cos 9 e | distribution depends strongly on the leptoquark mass. This is 
shown in Fig. ||]d where an additional distribution for a leptoquark mass of 80 GeV 
has been added. The cut |cos# e | < 0.8 is necessary to reduce the background 
from two-photon events, but it also significantly reduces the efficiency for small 
M LQ . 

Fig. |]a shows the selection efficiencies after all cuts as determined with ERATO- 
LQ as a function of the generated leptoquark mass Mlq. Fig. |]b shows the dis- 
tributions of Mjj for data, Standard Model background and for the state Vo using 
A = v^Tracm with M LQ = 80 GeV/c 2 and 120 GeV/c 2 . After all cuts, 21 events 
remain in the data which is in good agreement with the predicted 21.9±1.1 (stat) 
Standard Model background events. 

4.2 The neutrino plus hadronic jet channel 

In the case of the decay of a leptoquark into a neutrino and a single quark, the 
search has to be optimised for a single hadronic jet in the detector. Its transverse 
energy E T j et must be balanced by the neutrino. The cuts are therefore: 

• The event must contain more than four tracks (n c h > 4). 

• Exactly one jet must have been reconstructed in the pseudorapidity region 

tail < 2. 

• No hit in the MIP plug with a significant charge deposition is found. 

• The jet must consist of more than six particles (n^ > 6). 

• The distance of the primary vertex to the nominal interaction point has to 
be less than 2 cm. 

The following cuts are applied to the 432 data events which remain after this 
preselection: 

NQ1) The ratio between the jet energy E- ]et and the total visible energy E vis in the 
event has to be greater than 0.8 (Ej et /E vis > 0.8). This cut is very effective 
in reducing all sources of Standard Model background like multihadronic 
and two-photon events. 
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NQ2) The difference between the jet transverse energy, -Erjet, arid the missing 
transverse energy, has to be less than 2 GeV ( | -^Tjet — $r| < 2 GeV) 
in order to ensure that the final state consists of a single jet balanced by 
missing transverse energy. 



cuts 


Pre- 
selection 


(NQ1) 


Exjet — E T 

(NQ2) 


M LQ (GeV/c 2 ) 


state 


80 


scalar 


57.0% 


45.4% 


41.5% 


120 


62.3% 


54.4% 


48.3% 


160 


67.7% 


65.6% 


56.1% 


80 


vector 


59.0% 


43.6% 


39.9% 


120 


59.2% 


50.0% 


43.7% 


160 


62.2% 


60.1% 


51.6% 


e + e~ — > qq 


178.2 


2.9 


0.0 ±0.6 


e + e~ — > t + t~ 


4.6 


2.1 


0.3 ±0.1 


e + e~ — > 4 fermions 


71.0 


10.5 


6.3 ±0.3 


two-photon 


172.4 


4.6 


2.4 ± 1.1 


total BG 


426.2 


20.1 


8.9 ± 1.2 


data 


432 


24 


7 



Table 4: Selection efficiencies for three different leptoquark masses for scalar and 
vector leptoquarks. Also shown is the remaining number of data events and the 
expected number of background (BG) events after each cut of the neutrino plus 
hadronic jet selection. The Monte Carlo background is normalised to the data 
luminosity of 158.4 pb _1 at 189 GeV. The errors on the Monte Carlo background 
are statistical. 

Table |] shows the number of data and Monte Carlo events normalised to data 
luminosity after each cut beginning after the preselection. The signal efficien- 
cies for three different scalar and vector leptoquark masses are also given. The 
efficiencies are similar for scalar and for vector states. 

After all cuts, 7 events remain which is in good agreement with the Standard 
Model expectation of 8.9 ±1.2 (stat) events. From the two-photon events only the 
CC DIS events give a sizeable contribution to the final background composition 
together with the hadronic four-fermion processes. 

Fig. |5] shows some of the cut variables for data, Standard Model background 
and the leptoquark state S1/2 with a mass of 120 GeV/c 2 in arbitrary normali- 
sation. The sum of the Standard Model Monte Carlo distributions describes the 
data sufficiently well. Fig. |^a) shows the selection efficiencies as determined with 
the scalar and vector leptoquarks generated with ERATO-LQ after all cuts as a 
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function of the generated leptoquark mass Mlq. Fig. ||b) shows the distribution 
of the transverse mass Mp = after all cuts for data, Standard Model back- 
ground, and for the state S1/2 using A = i/47ra em with M L q = 80 GeV/c 2 and 
120 GeV/c 2 . 



5 Results 



The systematic errors on the expected signal rate are: (a) the luminosity mea- 
surement with less than 1 %, (b) the model dependence of the leptoquark frag- 
mentation with 1 to 5 %, (c) the parameter dependence for the jet finding with 
2 to 5 %, (d) the Monte Carlo statistics with 1 %. 

A special version of PYTHIA has been used to study the difference between 
models where the hadronisation is simulated before and after the leptoquark de- 
cay into an electron and a quark |38fl . The difference in the average charged 



multiplicity within the geometrical acceptance of the detector increases with lep- 
toquark mass and is always less than one unit. The model dependence of the 
leptoquark fragmentation has therefore been estimated by varying the cut on 
the charged multiplicity by one unit in the Monte Carlo while keeping it fixed 
in the data since the charged multiplicity cut is expected to be very sensitive 
to the hadronisation model. In addition, the jet finding parameters have been 
varied: the value of the minimum transverse energy for a jet has been changed 
by ±5 GeV and the jet radius R has been changed from 1 to 0.7. The effect from 
the variation of R is negligible. 

The Monte Carlo model dependence has been studied by comparing the al- 
ternative background Monte Carlo sets defined in section |3]. In the eq channel 
the alternative two-photon and four-fermion sample predict each an increase in 
the total number of background events of about 10%. This would lead to a 
higher limit on the leptoquark mass than the Monte Carlo sample used. For the 
i/q-channel the four-fermion generators yield consistent results and the expected 
two-photon background decreases from 2.4 to 1.3 events, leading to a negligible 
change in the mass limits. 

The systematic errors are added in quadrature and they are taken into account 
in the limit using the procedure of Highland and Cousins p9 |. 



The limits have been calculated for three different values of the branching 
ratio f3 = 1, P = 0.5, and (3 — > 0, since for (3 = no production in eq collisions is 
possible. The limit calculations have been performed according to the procedure 
of |I0] which takes into account the expected background, the mass resolution, the 
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candidates, and the efficiencies. The cross-section excluded at the 95% confidence 
level, (T95, resulting from these calculations is shown in Fig. [7| for scalar and 
for vector leptoquarks. The resulting mass limits are given in Table [5] for A = 
y/Ana cra , where the electromagnetic coupling constant a em is taken at the mass 
of the leptoquark with a em (MLQ) ~ 1/128. The upper limit at the 95 % CL of 
the coupling A (A') as a function of the mass Mlq is given in Figs. || and §. 



charge 


state 


(3=1 


(3 = 0.5 


(3 -> 


±1/3 


So, Si, S1/2, 


163 GeV/c 2 


158 GeV/c 2 


175 GeV/c 2 


-5/3 


Sl/2 


164 GeV/c 2 






-4/3 


So 


149 GeV/c 2 






-4/3 


Si 


156 GeV/c 2 






-2/3 


Sl/2) S1/2, S1/2 


121 GeV/c 2 




141 GeV/c 2 


+2/3 


Si 






162 GeV/c 2 



Table 5: Mass limits for scalar leptoquarks and squarks for A = ^Aira^ and the 
different (3 values. 



We have also derived the limit on Mlq as a function of (3 for couplings of 
electromagnetic strength (A = ^4ira cm ), i.e. the assumption is dropped that the 
branching ratio (3 can take only the values 0,0.5 and 1. In Fig. [12| the limit on 



(3 as a function of M L q is compared to the DO results |L4| which only exclude 
values up to M L q = 80 GeV/c 2 for (3 = 0. Our analysis is more sensitive in the 
low (3 region, yielding a limit of Mlq = 175 GeV/c 2 for (3 — > 0. A similar region 
has recently been explored by HI 



ro 



Exactly the same procedure as for scalar leptoquarks has been used to deter- 
mine the mass limits and limits on the Yukawa coupling for vector leptoquarks. 
The systematic errors are also identical. The results for the mass limits are shown 
in Table 0. The upper limit at the 95 % CL of the coupling A as a function of 



the mass Mlq is given in Figs. 1 1 0|— (13 



5.1 Conclusions 



We have searched for singly-produced leptoquarks in electron-photon interactions 
at an e + e~ centre-of-mass energy of 189 GeV using the data collected with the 
OPAL detector at LEP. The data correspond to an integrated luminosity of about 
160 pb -1 . Some scalar leptoquark states can also be identified with squarks in 
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charge 


state 


= 1 


/3 = 0.5 





-1/3 


Vl/2, Vx/2 


176 GeV/c 2 




182 GeV/c 2 


+ 1/3 


Vi 






188 GeV/c 2 


-5/3 


V 


177 GeV/c 2 






-5/3 


Vi 


182 GeV/c 2 






-4/3 


Vi /2 


152 GeV/c 2 






±2/3 


V , Vx, V 1/2 


151 GeV/c 2 


149 GeV/c 2 


163 GeV/c 2 



Table 6: Mass limits for vector leptoquarks for A = ^4na cm and the different 
possible (5 values. 



R-parity violating SUSY models. No evidence was found for the production of 
these particles. Therefore, 95% confidence limits were set on the mass as well as 
on the Yukawa coupling A for scalar and vector leptoquarks and A' for squarks 
as a function of the mass for different branching fractions f3 in eq final states. 



Acknowledgements 



We thank M. Doncheski, S. Godfrey and U. Katz for very useful discussions and 
C. Papadopoulos for his help in using ERATO-LQ. 

We particularly wish to thank the SL Division for the efficient operation of the 

LEP accelerator at all energies and for their continuing close cooperation with 

our experimental group. We thank our colleagues from CEA, DAPNIA/SPP, CE- 

Saclay for their efforts over the years on the time-of-flight and trigger systems 

which we continue to use. In addition to the support staff at our own institutions 

we are pleased to acknowledge the 

Department of Energy, USA, 

National Science Foundation, USA, 

Particle Physics and Astronomy Research Council, UK, 

Natural Sciences and Engineering Research Council, Canada, 

Israel Science Foundation, administered by the Israel Academy of Science and 

Humanities, 

Minerva Gesellschaft, 

Benoziyo Center for High Energy Physics, 

Japanese Ministry of Education, Science and Culture (the Monbusho) and a grant 
under the Monbusho International Science Research Program, 
Japanese Society for the Promotion of Science (JSPS), 
German Israeli Bi-national Science Foundation (GIF), 



17 



Bundesministerium fur Bildung und Forschung, Germany, 
National Research Council of Canada, 
Research Corporation, USA, 

Hungarian Foundation for Scientific Research, OTKA T-029328, T023793 and 
OTKA F-023259. 



References 

[1] W. Buchmuller, R. Riickl and D. Wyler, Phys. Lett. B191 (1987) 442. 

[2] J. Kalinowski, R. Riickl, H. Spiesberger and P.M. Zerwas, Z. Phys. C74 
(1997) 595. 

[3] H.P Nilles, Phys. Rep. 110 (1984) 1; 

H.E. Haber and G.L. Kane, Phys. Rep. 117 (1985) 75. 

[4] M.A. Doncheski and S. Godfrey, Phys. Lett. B393 (1997) 355; 

M.A. Doncheski and S. Godfrey, Mod. Phys. L. A12 (1997) 1719 and private 
communications. 

[5] T.M. Aliev, E. Iltan and N.K. Pak, Phys. Rev. D54 (1996) 4263. 

[6] G. Belanger, D. London and H. Nadeau, Phys. Rev. D49 (1994) 3140. 

[7] J. Bliimlein and E. Boos, Nucl. Phys. Proc. Suppl. B37 (1994) 181. 

[8] J.L. Hewett, S. Pakvasa, Phys. Lett. B227 (1987) 178. 

[9] DELPHI Collab., P.P. Allport et al, Phys. Lett. B446 (1999) 62. 

[10] HI Collab., C. Adloff et al, Eur. Phys. J. Cll (1999) 447. 

[11] ALEPH Collab., R. Barate et al, Eur. Phys. J. C12 (2000) 183; 
OPAL Collab., G. Abbiendi et al, Eur. Phys. J. C6 (1999) 1; 
L3 Collab, M. Acciari et al, Phys. Lett. B433 (1998) 163. 

[12] DELPHI Collab., P. Abreu et al., Phys. Lett. B316 (1993) 620; 
DELPHI Collab., P. Abreu et al., Phys. Lett. B275 (1992) 222; 
ALEPH Collab., D. Decamp et al., Phys. Rept. 216 (1992) 253; 
OPAL Collab., G. Alexander et al., Phys. Lett. B263 (1991) 123; 
L3 Collab, B. Adeva et al., Phys. Lett. B261 (1991) 169. 

[13] OPAL Collab., G. Abbiendi et al, Eur. Phys. J. C13 (2000) 15. 



18 



[14] DO Collab., B. Abbott et al, Phys. Rev. Lett. 83 (1999) 2896; 
DO Collab., B. Abbott et al., Phys. Rev. Lett. 80 (1998) 2051. 

[15] CDF Collab., F. Abe et al, Phys. Rev. Lett. 82 (1999) 3206; 

CDF Collab., F. Abe et al., Phys. Rev. Lett. 81 (1998) 4806; 

CDF Collab., F. Abe et al, Phys. Rev. Lett. 79 (1997) 4327; 

CDF Collab., F. Abe et al, Phys. Rev. Lett. 75 (1995) 1012 

[16] C. Grosso-Pilcher et al., Combined Limits on First Generation Lepto- 
quarksfrom the CDF and DO Experiments, FERMILAB-PUB-98-312-E, [Eep] 
| ex/98100T5| . 

[17] DO Collab., B. Abbott et al., Phys. Rev. Lett. 83 (1999) 4476; 
DO Collab., F. Abe et al, Phys. Rev. Lett. 83 (1999) 2133; 
CDF Collab., B. Abbott et al., Phys. Rev. Lett. 83 (1999) 4937; 
CDF Collab., F. Abe et al., Phys. Rev. D56 (1997) 1357. 

[18] ALEPH Collab., R. Barate et al., Search for R-Parity Violating Decays 
of Supersymmetric Particles in e + e~ Collisions at Center- of- Mass Energies 
from 189 GeV to 202 GeV, CERN-EP-2000-132; 
DELPHI Collab., P. Abreu et al., Eur. Phys. J. C13 (2000) 591; 
OPAL Collab., G. Abbiendi et al, Eur. Phys. J. C12 (2000) 1; 
L3 Collab., M. Acciarri et al., Phys. Lett. B459 (1999) 354. 

[19] OPAL Collab., K. Ahmet et al., Nucl. Instr. and Meth. A305 (1991) 275. 

[20] S. Anderson et al, Nucl. Instr. and Meth. A403 (1998) 326. 

[21] OPAL Collab., K. Ackerstaff et al., Eur. Phys. J. C12 (2000) 551. 

[22] B.E. Anderson et al., IEEE Transactions on Nuclear Science 41 (1994) 845. 

[23] OPAL Collab., K. Ackerstaff et al., Phys. Lett. B391 (1997) 221. 

[24] P. Aurenche and G. A. Schuler (conveners), 77 Physics, in "Physics at 
LEP2", CERN 96-01, eds. G. Altarelli, T. Sjostrand and F. Zwirner, Vol. 1 
(1996) 291. 

[25] C. Papadopoulos, Comp. Phys. Comm. 118 (1999) 118; 
C. Papadopoulos, Comp. Phys. Comm. 101 (1997) 183. 

[26] T. Sjostrand, Comp. Phys. Comm. 82 (1994) 74. 

[27] M. Gliick, E. Reya and A. Vogt, Phys. Rev. D46 (1992) 1973; 
M. Gliick, E. Reya and A. Vogt, Phys. Rev. D45 (1992) 3986. 

[28] G.A. Schuler and T. Sjostrand, Z. Phys. C68 (1995) 607. 



19 



[29] S. Jadach, B.F.L. Ward and Z. Was, Comp. Phys. Comm. 79 (1994) 503. 

[30] S. Jadach, W. Placzek, B.F.L. Ward, in 'Physics at LEP2', eds. G. Altarelli, 
T. Sjostrand and F. Zwirner, CERN 96-01, vol.2 (1996). 

[31] G. Marchesini, B.R. Webber, G. Abbiendi, I.G. Knowles, M.H. Seymour and 
L. Stanco, Comp. Phys. Comm. 67 (1992) 465. 

[32] R. Engel and J. Ranft, Phys. Rev. D54 (1996) 4244; 
R. Engel, Z. Phys. C66 (1995) 203. 

[33] J.A.M. Vermaseren, Nucl. Phys. B229 (1983) 347. 

[34] A. Buijs, W.G.J. Langeveld, M.H. Lehto and D.J. Miller, Comp. Phys. 
Comm. 79 (1994) 523. 

[35] J. Fujimoto et al., in "Physics at LEP2", eds. G. Altarelli, T. Sjostrand and 
F. Zwirner, CERN 96-01, Vol. 2 (1996). 

[36] S. Jadach, W. Placzek, M. Skrzypek, B.F. Ward and Z. Was, Corn- 
put. Phys. Commun. 119 (1999) 272. 

[37] OPAL Collab., K. Ackerstaff et al., Z. Phys. C73 (1997) 433. 

[38] C. Friberg, E. Norrbin, T. Sjostrand, Phys. Lett. B403 (1997) 329 and private 
communications. 

[39] R.D. Cousins and V.L. Highland, Nucl. Inst, and Meth. A320 (1992) 331. 

[40] P. Bock, Determination of exclusion limits for particle production using dif- 
ferent decays channels with different efficiencies, mass resolutions and back- 
grounds, Heidelberg University, HD-PY-96/05 (1996). 



20 




q 



Figure 1: Diagram of the s-channel production of a leptoquark (LQ) in electron- 
photon scattering. The photon is radiated by one of the LEP beams, fluctuates 
into a hadronic object and one of the quarks interacts with an electron from the 
other beam. 




(b) 

Figure 2: Alternative single leptoquark production mechanisms, where the pho- 
ton interacts pointlike: (a) "direct" interaction with the beam electron, (b) pho- 
ton is absorbed by the leptoquark "emitted" by the other beam electron. 
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Figure 3: Electron-quark decay channel: (a) distribution of the ratio (E^ + 
E e )/E vis after the preselection; (b) distribution of the ratio /M^ after ap- 
plying the additional cut EQ1; (c) distribution of the difference \E e — E e ^\ after 
applying the additional cut EQ2; (d) distribution of the electron scattering angle 
| cos# c | after applying cuts EQ1 through EQ3 to the preselected data. The points 
with error bars are the data, the full line represents the total Standard Model 
background normalised to data luminosity and the dashed (dotted) histogram 
shows the distribution for the vector leptoquark state Vo with a mass of 120 
GeV/c 2 (80 GeV/c 2 ). The normalisation of the leptoquark signals is arbitrary. 
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Figure 4: (a) Selection efficiency in the eq channel after all cuts for scalar (dots) 
and vector (squares) leptoquarks; (b) Invariant mass Mjj of the two jets (= 
leptoquark mass) after all cuts for data (points with error bars), Standard Model 
background (full line) and two different vector leptoquark masses (dotted line 
for 80 GeV/c 2 and dashed line for 120 GeV/c 2 ). The state Vo was chosen to 
normalise the signal, using A = y/Ana^. 
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Figure 5: Neutrino-quark decay channel: (a) distribution of the ratio Ej et /E v i s 
after the preselection; (b) distribution of the variable |Ej et — $r | after the cut 
NQ1 ((Ej et /E vis > 0.8). The points with error bars are the data, the full line 
represents the total Standard Model background normalised to data luminosity 
and the dashed histogram shows the distribution for the scalar state S1/2 with a 
mass of 120 GeV/c 2 . The normalisation of the leptoquark signals is arbitrary. 

23 




Figure 6: (a) Selection efficiency in the vq channel after all cuts for scalar (dots) 
and vector (squares) leptoquarks; (b) transverse mass Mp — after all cuts 
for data (points with error bars), Standard Model background (full line) and two 
different scalar leptoquark masses (dotted line for 80 GeV/c 2 and dashed line for 
120 GeV/c 2 ). The state Si/ 2 (charge -2/3) was chosen to normalise the signal, 
using A = y/4na cm . 



J2 
a. 


0.45 


ON 


0.4 


t> 






0.35 




0.3 




0.25 




0.2 




0.15 




0.1 




0.05 








scalar leptoquarks 



1 1 1 1 1 1 1 1 1 11 1 


| ' ' ' ' | ' ': 


r(a) OPAL \ 




r — eq channel \ 




r ■ - - vq channel \ 


S 1/2 (-5/3) " 


T both channels 




j^j 

- 1 1 1 1 1 1 1 1 1 1 1 1 1 


" "'A. ^Si: 

1 i i i i 1 \_- 



100 125 150 175 



■|0.45 

8 0.4 

b 

0.35 
0.3 

0.25 
0.2 

0.15 
0.1 

0.05 




vector leptoquarks 



I 1 1 1 1 I 1 1 
"(b) OPAL 

— eq channel 

■ - - vq channel 

both channels 




M LQ [GeV/c] 



M LQ [GeV/c 2 ] 



Figure 7: Cross-section excluded at the 95% confidence level, 095, using the 
number of candidates in the data for each channel separately, corresponding to 
(3=1 for the electron-quark channel (full line) and (3=0 for the neutrino-quark 
channel (dashed line), as well as for equal branching ratio into both channels 
((3=0.5, dotted line). The expected SM background, the mass resolution, the 
candidates, and the efficiencies are taken into account in the calculation. The 
cross-section for the production of the states a) Si/ 2 (charge —5/3) and b) V are 
also shown. 
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Figure 8: Limits on the coupling constants A (A' ljfc ) for scalar leptoquark states 
(squarks) with (a) charge ±|: S with (3 — 1 (full line), S , <1r, Si with (3 = 0.5 
(dashed line) as well as Si/ 2 and d L with (3 — (b) charge — |: Si/ 2 with (3 = 1. 
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Figure 9: Limits on the coupling constant A (A^) for scalar leptoquark states 
(squarks) with (a) charge — |: Si with (3 — 1 (full line), S , with f3 = 1 (dashed 
line) (b) charge ±|: Si/ 2 , u L , Si/ 2 with (3 = 1 (full line), Si/ 2 with (3 = (dotted 
line) and Si with (3 = (dashed-dotted line). 
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Figure 10: Limits on the coupling constant A for vector leptoquark states with 
(a) charge ±|: Vi /2 ,Vi /2 with (3 = 1 (full line) and V 1/2 , V x with (3 = (b) 
charge — |: V with (3 = 1 (full line) and Vi with (3 = 1 (dashed line). 
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(a) OPAL, Vs ee =189 GeV 
charge -4/3, 95% CL excl. 
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Figure 11: Limits on the coupling constant A for vector leptoquark states with 
(a) charge -§: V 1/2 with (3 = 1 (full line) (b) charge ±|: V with (3 = 1 (full 
line), V , Vi with (3 = 0.5 (dashed line) and Vi/ 2 with (3 = (dotted line). 
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Figure 12: Limit on Mlq as a function of (3 for couplings of electromagnetic 
strength (A = y/^Q. CTD ) for the charge ±1/3 scalar states. The limit on /3 is 



compared to the DO results |p^ . For (3 = no production in eq collisions is 
possible. 
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